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Abstract 

The intrinsic Zeeman energy is precisely one half of the cyclotron energy for electrons in graphene. As 
a result a Landau-level mixing occurs to create the energy spectrum comprised of the 4j-fold degenerated 
zero-energy level and 4-fold degenerated nonzero-energy levels in the j-layer graphene, where j — 1,2,3 for 
monolayer, bilayer and trilayer, respectively. The degeneracy manifests itself in the quantum Hall (QH) 
effect. We study how the degeneracy is removed by the Coulomb interactions. With respect to the zero- 
energy level, an excitonic gap opens by making a BCS-type condensation of electron-hole pairs at the filling 
factor V = Q. It gives birth to the Ising QH ferromagnet aXv — ±1 for monolayer, v = ±1, ±3 for bilayer, and 
V = ±1,±3,±5 for trilayer graphene from the zero-energy degeneracy. With respect to the nonzero-energy 
level, a remarkable consequence is derived that the effective Coulomb potential depends on spins, since 
a single energy level contains up-spin and down-spin electrons belonging to different Landau levels. The 
spin-dependent Coulomb interaction leads to the valley polarization at v = ±4, ±8, ±12, • • • for monolayer, 
ly — ±2, ±6, ±10, • • • for bilayer, and ly = ±2, ±4, ±8, ±12, • • • for trilayer graphene. 
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I. INTRODUCTION 



Recent experiments has established that the charge carriers in graphene are massless Dirac 
electronsii^i^i^i^. Dirac electrons have a peculiar property that they have the intrinsic Zeeman energy 
precisely one half of the cyclotron energy in magnetic field. Consequently a Landau level mixing 
occurs so that one energy level contains up-spin and down-spin electrons coming from different 
Landau levels. It has two important consequences; the emergence of the zero-energy stated, and 
the degeneracy of the up-spin and down-spin states for each nonzero-energy level [FiglUa)]. The 
aim of this paper is to explore new phenomena due to this intrinsic Zeeman effect in graphene. 

The quantum Hall effect (QHE) in graphene^'^'^'^i^ is unconventional. The filling factors form a 
series [FiglUb)], 

u = 0, ±1, ±2, ±4, ±6, ±8, ±10, ±12, • • • , (1.1) 

where the bold-face series had been predicted^'^'^ before it was found experimentallyii^i^i^, while 
the full series was discovered latep^ with larger magnetic field applied. Subsequently theoretical 
works^'ii'^'^'^ have been made to interpret the series. Furthermore, the series reads 

1/ = 0, ±1, ±2, ±3, ±4, ±6, ±8, ±10, ±12, • • • , (1.2) 

in a bilayer graphene, where the bold-face series has been found experimentally^ and studied 
theoretically^^, while it reads 

1/ = 0, ±1, ±2, ±3, ±4, ±5, ±6, ±8, ±10, ±12, • • • , (1.3) 

in a trilayer graphene, where the bold-face series has been predicted theoretically^. In this paper 
we show that the full series emerge when Coulomb interactions become important. It is notable 
that the basic height in the Hall conductance is 4e^//i for the bold-face series for all these graphene 
systems, indicating the 4-fold degeneracy of the energy level except for the first step at the z/ = 
point within noninteracting theory. 

Conduction and valence bands in graphene form conically shaped valleys, touching at a point 
[Figl2]. There are two inequivalent Brillouin zone corners, called the K and K' points, at which 
massless Dirac electrons emerge^'^'^. Let us refer to the valley as the Dirac valley, and assign the 
valley index to the electron so that the electron at the K (K') point carries the index r = ± (— ). As 
we have mentioned, the graphene QH system is characterized by the emergence of the zero-energy 
state and the degeneracy of the up-spin and down-spin states for each nonzero-energy level. Since 
this holds separately at the K and K' points, each energy level has a 4-fold degeneracy [Fig|l](a)], 
and the noninteracting system has the SU(4) symmetry. 

Intriguing phenomena occur when we introduce Coulomb interactions. It is necessary to consider 
the Coulomb effect in the zero-energy state and nonzero-energy states separately. 

The zero-energy state is distinctive, since it contains both electrons and holes. Electron-hole 
pairs form an excitonic condensation due to attractive interaction, producing an excitonic gap. We 
obtain a BCS-type state of electron-hole pairs at u = 0. As a consequence of excitonic condensation, 
the degeneracy of the zero-energy state is resolved into two subbands each of which contains either 
electrons or holes. The Coulomb Hamiltonian, projected to each of them, has the U(l) symmetry 
but is broken into the Z2 symmetry. We obtain the Ising QH ferromagnet at u = ±1. 

We next study nonzero-energy states of electrons, where a single energy level contains up-spin 
and down-spin electrons belonging to different Landau levels [FiglUa)]. We derive a remarkable 
consequence that the effective Coulomb potential depends on the spin and the valley in each energy 
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FIG. 1: (Color online) (a) The energy level (number in circle) and the Landau level (number with index E 
or H) for electrons and holes. The spin is indicated by a solid red (open blue) arrow for electron and hole 
at the K and K' points. The iVth energy level contains up(down)-spin electrons from the A^th Landau level 
and down(up)-spin electrons from the (7V-l)th Landau level at the K (K') point for > 1. The zero-energy 
level {N — 0) contains up(down)-spin electrons and down(up)-spin holes at the K (K') point, (b) The QH 
conductivity in graphene. The dotted black curve shows the sequence v — ±2, ±6, ±10, • • • , while the solid 
red curve the sequence ^ = 0, ±1, ±2, ±4, • • • . 



leveL Namely, it depends on the spin of electrons as well as on which valley electrons belongs to. 
The Coulomb Hamiltonian, projected to a single energy level, possesses only the U(1)®U(1)(8'Z2 
symmetry. Consequently, the Coulomb interaction resolves the 4-fold degeneracy into two 2-fold 
degeneracies, generating a new series at = ±4, ±8, ±12, • • • . The ground state is a valley polarized 
state, by which we mean that more electrons are present in one valley than the other. This occurs 
because the Coulomb energy is lower in higher Landau levels. 

Our analysis can be generalized to the j-layer graphene system, where j = 1, 2 and 3 correspond 
to monolayer, bilayer and trilayer, respectively. We show that the zero-energy states have the 
4j-fold degeneracy. As a result the quantized values of the Hall conductivity become 

fTxy = ±(n±0^, n = 0,1,2,- ••, (1.4) 

within noninteracting theory. We also discuss the effects due to Coulomb interactions. 

This paper is composed as follows. In Section [U we investigate the energy spectrum in the 
presence of magnetic field. We show that two Landau levels mix to create one energy level. In 
Section IIIII we construct the projected Coulomb Hamiltonian relevant to analyze physics in the 
A'^th energy level. The effective Coulomb potential is shown to depend on the spin and the valley 
through the form factors characterizing Landau levels. In Section IIVI we treat the zero-energy 
level {N = 0), which contains both electrons and holes. We present a clear-cut description of an 
excitonic condensation, producing a gap to the electron and hole states, on the analogy of the BCS 
superconductor. In Section |Vl studying how Coulomb interactions resolve the degeneracy of the 
nonzero-energy level, we derive the valley polarized ground state. In Section [Vll our analysis is 
generalized to multilayer graphene systems. Finally, Section IVIIl is devoted to discussions. 
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FIG. 2: (Color online) (a) Graphene is made of a honeycomb lattice. It consists of two triangular sublattices 
generated by the two basis vectors ai and 02 from the base points A and B in primitive cell (dotted rectangle), 
(b) The reciplocal lattice is also a honeycomb lattice. The first Brillouin zone is depicted together with the 
Dirac valleys describing the low-energy band structure. 

II. DIRAC ELECTRONS AND INTRINSIC ZEEMAN EFFECT 
A. Dirac Hamiltonian 

The low-energy band structure of graphene is described by conically shaped valleys, touching at 
a point. There exists two inequivalent Brillouin zone corners ±K = (±47r/\/3a, 0) with a the lattice 
constant. They are called the K and K' points [Figl2lb)]. The two-dimensional energy dispersion 
relation is linear in these Dirac valleys, 

£r{k) = hVF\k-TK\, (2.1) 

where vp is the Fermi velocity, and r = it is the valley index. This linear behavior has been 
confirmed experiment ally^S up to 3eV. It is convenient to make the change of variable, k = k' + tK, 
and rewrite the dispersion relation as 

S{k') = hvp\k'\. (2.2) 

We use the variable k' to show that hk' is the relative momentum of electrons, |fc'| <C \K\, measured 
from tK at the K or K' point. 

There exists one electron per one carbon and the band-filling factor is 1/2 in graphene. Hence 
we have electron excitations in the conduction band and hole excitations in the valence band. 
Furthermore there exists the electron-hole symmetry. 

The dispersion relation (12.2(1 is that of 'relativistic' Dirac fermion a^^'^^'^^ . Hence the second- 
quantized Hamiltonian is given by 

H =^ I (fx-^l{x)Hl^^r{x), (2.3) 
where is the quantum mechanical Hamiltonian, 

= {tOxPx + CfyPy) 75> (2-4) 

with Pi = —ihdi, the Pauli matrix Ui and the Dirac matrix 75. The field ^ri'x) consists of eight 
components, corresponding to the spin degree of freedom, the electron-hole degree of freedom and 
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the valley degree of freedom. Being massless, the Hamiltonian describes Weyl fermions, and the 
spin index a represents the helicity. 

The field operator '^t{x) is expanded in terms of the eigenfunctions of as 

^r{x) = ^er{x) + ^-hrCa;), (2.5) 

with 

*e.(^) = E / |^<(fc')<(fc')e^'"'^ (2.6a) 

M'h.(x) = E / (2.6b) 

where c%(h!') and dp^k') are the annihilation operator of an electron and the creation operator 
of a hole, respectively, while u'^{k') and v'^{k') are the corresponding eigenfunctions of the Dirac 
operator (|2.4p with the eigenvalue S{k') = ±ihv-p\k'\. 

The field ^t{x) describes solely the dynamics of electrons localized in each Dirac valley at the 
K or K' point. In the dispersion relation (|2.2p . Kk' is the momentum fluctuation around zizhK 
at the K or K' point. The total momentum of the electron is not hk' but hk = hk' it hK. Due 
to this fact, the eigenfunctions u^ik') and v^{k') are not wave functions but envelope functions^. 
The wave functions are given by e^'^^''^u'^{x) and e^'^^''^v'^{x). Accordingly, the field operators of 
electrons and holes are 

tl;,Ax) = e*-^-^M/e.(^r), V'hr(a^) = e'"^-"^hr(a;). (2.7) 

Thus, to analyze the Coulomb interaction, it is necessary to use ■ip^^{x) and ■i/'hrC^)- Appendix 
Owe derive the Dirac Hamiltonian (12. 4|) based on the effective-mass description^. 

We make a comment on the factor r in front of ax in (12. 4p . It has the standard expression of the 
Dirac Hamiltonian for r = +, that is, at the K point. The K point is transformed into the K' point 
under the mirror reflection. Corresponding to this we have the mirror-reflected Dirac Hamiltonian 
(|2.4p for r = — , that is, at the K' point. The total Hamiltonian (|2.3I) has the mirror symmetry. 



B. Landau Levels 

We apply the magnetic field to a graphene sheet taken on the xy plane. It is introduced to the 
Hamiltonian (|2.3|) by making the minimal substitution, 

= VF {ra^Px + (TyPy) Ts, (2-8) 

where Pi = —ihdi + eAi is the covariant momentum. We assume a homogeneous magnetic field 
S = V X A = (0, 0, —B) with B > along the z axis. The presence of the external magnetic 
field modifies the mirror symmetry as follows: The modified mirror reflection not only transforms 
the K point into the K' point but also reverses the direction of the magnetic field to maintain the 
symmetry. This is most clearly seen in ()2.19|) . as we shall discuss later. Consequently, the K and 
K' points become physically distinguishable. 

Electrons make cyclotron motion in magnetic field. The helicity is no longer a good variable, 
since there exists a special direction for the spin, that is, the direction of magnetic field. In this 
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case it is convenient to use the standard representation for the Dirac matrices, where 



75 



1 

1 



(2.9) 



In order to discuss the QHE, we solve the quantum mechanical problem with the Hamiltonian (|2.8I 
in each Dirac valley, 



and we expand the field operator in terms of the new eigenfunctions. 
To solve the eigen equation (|2.10p we express the Hamiltonian (|2.8|) as 



(2.10) 



Qr 
Qr 



with the use of (12. 9p for 75, where 



It is diagonalized as 



ifp = diag. ( \/QTQT,-\/QTQr] ■ 



We introduce a pair of operators 



iBiPcc + iPy) 



^{Px-iPy 



V2h 



V2h 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



satisfying [o, a^] = 1, where £3 = \/h/eB is the magnetic length. Since the operators (I2.12P is 
rewritten as 



at \ a 

a \ at 



(2.15) 



with LUc = V^fivp/^B, the diagonalized Hamiltonian reads 



Ht^ =hjjja diag. ( V ata, V aat, — v ata, —V aat 



=hwc diag 



aat, V Qrta — vaa',— vaial. 



(2.16a) 
(2.16b) 



Since ata is the number operator, the eigenvalue of the Dirac Hamiltonian follows immediately. 



£+ = [Vn, VnTT, - Viv, -VnTT) , 



(2.17a) 
(2.17b) 
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with the eigenstate 



\N) = -^{aYm, (2.1^ 



where A^ = 0,1,2, 3,---. 

The operators a and are the Landau-level ladder operators. Hence, N in (|2.17p represents 
the Landau-level index. We have found that the energy of an electron in the A^th Landau level is 
either HlOcVN or fiWc\/N + 1. It is curious that the energy of an electron in the lowest Landau 
level {N = 0) is zero though it performs cyclotron motion. This puzzle is solved in the succeeding 
subsection. 



C. Pauli Hamiltonian 

To reveal the intrinsic structure of the energy spectrum, we investigate the Hamiltonian- 



= Q±Q± = vl {-ihV + eAf =feha^B , (2.19) 

which is the building blocks of the Dirac Hamiltonian (|2.13p . Here, the direction of the magnetic 
field is found to be effectively opposite at the K and K' points. Since this has the same form as 
the Pauli Hamiltonian with the mass m* = l/4fp except for the dimension, we call it the Pauli 
Hamiltonian for brevity. The salient feature of the relativistic Dirac Hamiltonian is that its spectrum 
is mapped from that of the nonrelativistic Pauli Hamiltonian. Thus, the energy eigenvalue £J^ of 
the Dirac Hamiltonian is constructed as = zty^E^ from the energy eigenvalue EJ^ of the 
Pauli Hamiltonian Hp. 

In the Pauli Hamiltonian (12.191) . the first term is the kinetic term while the second term is the 
Zeeman term. It is fixed uniquely as an intrinsic property of the Dirac theory: We call it the 
intrinsic Zeeman effect. 

The Landau level is created by electrons making cyclotron motion. In the conventional QHE, 
since the Zeeman energy can be considered much smaller than the Landau-level separation, we 
may treat it as a perturbation. However, this is not the case in graphene. According to the Pauli 
Hamiltonian (|2.19|) . the intrinsic Zeeman energy is precisely one half of the cyclotron energy for 
Dirac electrons, and two Landau levels mix to create one energy level, as illustrated in Fig|l](a). 

We consider the K point (r = -|-). It is obvious that the up-spin and down-spin states are 
eigenstates of the Pauli Hamiltonian (I2.19P and hence eigenstates of the Dirac Hamiltonian (I2.13p . 
and that the up-spin state has a lower energy than the down-spin state when they belong to the 
same Landau level. On the other hand, the direction of the spin is opposite at the K and K' points. 
Hence, we can make the following identification of quantum numbers for electrons in the A^th level 

> 1) of the energy spectrum [FiglUa)], 

= 4ii = ^N^ = ^nU = f^c^- (2.20) 

A similar identification can be made for holes in the A^th energy level {N > 1) . The zeroth energy 
level (A^ = 0) consists of the up-spin electron and the down-spin hole at the K point, and the 
down-spin electron and the up-spin hole at the K' point [FiglUa)], coming from the lowest Landau 
levels for electrons and holes. 

Counting the states at the K and K' points all together, one energy level has a 4-fold degener- 
acy. Each filled energy level contributes one conductance quantum e'^/h to the Hall conductivity. 
Consequently the resulting series is = ±2, ±6, ±10, • • • , which accounts for the bold-face series 
(|l.ip in the monolayer graphene. 
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D. Field Operators 



We focus on electrons in the A''th energy level. (Essentially the same analysis is applicable to 
holes.) We start with the quantum-mechanical state in a single Landau level. We decompose the 
electron coordinate x = {x,y) into the guiding center X = {X,Y) and the relative coordinate 
R = {R^,Ry), X = X + R, where R^ = -Py/eB and Ry = Px/eB with P = {Px,Py) the 
covariant momentum. They satisfy the commutation relations [X, ^] = —1^%, [Px,Py] = i^^/^'s, 
[X, Px] = [X, Py] = [y, Px] = [y, Py] = 0. We define a set of operators 6, 6^ by 

obeying \b, 6^] = 1, in addition to a set of operators a, a) by (I2.14p . The quantum-mechanical states 
are the Fock states, 

|iV,n) = ^i=(Gt)^(&t)"|o) (2.22) 
V A* !n! 

with |0) the Fock vacuum, a|0) = 6|0) = 0. 

We have constructed the energy spectrum of the Dirac Hamiltonian (|2.13l) in each Dirac valley, 
as illustrated in FiglUa). The A^th energy level {N ^ 0) contains electrons coming from two Dirac 
valleys (r = it), whose field operators are expanded in terms of the eigenfunctions. Corresponding 
to "^Tv-i' ^N^ "^TV-i ™ (|2.20l) . we have 

<+(x) = e^^" Y.^x\N, n)c\{N, n), (2.23a) 

n 

i^N-{x) = e-'^'^^{x\N,n)cl{N,n), (2.23b) 

n 

i^N+i.x) = e*-^^ ^{x\N - 1, n)c^(7V - 1, n), (2.23c) 

n 

^\_{x) = e~^'^^^(a;|iV- l,n)cl(iV-l,n), (2.23d) 

n 

where c%{N,n) is the annihilation operator acting on the Fock state \N,n) at the r point. In what 
follows we suppress the Landau-level index in c%{N,n) with the understanding that c\{n) = 
c\{N,n), ci(n) = ci(iV,n), c^(n) = c:'^(iV - 1, n), cl(n) = cI(iV- l,n). See with respect 

to the factor e^^^'^ . The electron operators in the zeroth energy level are given by ■i/'jv+(^) ^^^d 
i^N-ix) with = 0. 



III. COULOMB INTERACTIONS 



We make the basic assumption that the cyclotron energy is much larger than the Coulomb energy 
as in the conventional QHE^^i^^i^I. We neglect the Landau-level mixing by the Coulomb interaction, 
and analyze it within one energy level [FiglUa)]. There exists a consistent formalism, known as the 
Landau-level projection. The projected theory presents not only a good approximation but also 
an essential way to reveal a new physics inherent to the QH system24i2^i2I. Then we may treat 
the electron system, the hole system and the zero-energy system independently of each other. The 
hole system has the same structure as the electron system due to the electron-hole symmetry. On 
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the other hand, the zero-energy system contains both electrons and holes. We study the electron 
system in this section, and the zero-energy system of electrons and holes in Section HVl 

Assuming all lower levels are filled up, we study the Coulomb interaction between electrons 
confined within the A'^th energy level {N > 1). The Coulomb Hamiltonian reads 

H=\j j d'xd^y V{x - y)p{x)p{y), (3.24) 

where p {x) is the density operator 

pix) = ^rrH^)rr'i^)- (3.25) 
rT'(T 

The Landau-level projection of the Coulomb Hamiltonian (I3.24p is a simple generalization of the 
lowest-Landau-level projections^ familiar in the conventional QHE. We require the density operator 
to be comprised solely of the electron fields belonging to the A^th energy level. Thus, from (|3.25p 
we construct the projected density operator Pi\f{x) by 

Pn{^) = Yl rJri^WNr'i^), (3-26) 
tt'ct 

where 'i/'Arr(^) the field operator (|2.23ll . As we show in Appendix HTl the projected density 
operator is rewritten in the momentum space as 



Pn 



[q) = Y,F:^^,{q)bZ{q). (3.27) 



where D'^'^i (q) is the bare density operator—, 

D^f,{q) = ^ y (m|e-*[''+^-^-"'^l^|n)<t(^)c-:(n), (3.28) 



27r 

mn 



with 



and F^^t (g) is the form factor, 

(3.29a) 
(3.29b) 
(3.29c) 
(3.29d) 

(3.30a) 
(3.30b) 

Here we have used the relation, q = q^ 31^, due to the lattice structure. It should be remarked 
that the bare density operator D'^!^, (q) involves only the guiding center coordinate X, while the 
form factor F^^, (q) involves only the relative coordinate R. 
The form factors are explicitly given by using^ 

V(iV + Af)! VV2; " \ 2 J 



4+ i1) 


= F^_ 


(q) = Fn (q) , 




= fI_ 


(q) = Fn-1 (q) 




= Fi+ 


(q) = Gn (q) , 


(q) 


= fI_ 


(q) = i-q) 


FNiq) 


={N\e~ 




Gn (q) 


= {N\e- 
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FIG. 3: (Color online) The spin dependence of the effective Coulomb potential. The vertical axis is the 
potential in unit of e'^/AneiB- The horizontal axis is the distance r in unit oi £b- The dotted red curve 
represents the ordinary Coulomb potential without the form factor, (a) The effective Coulomb potential 
Vl\r) in the iVth energy level for N = 0,1,2,3 from top to bottom. Note that V^^{r) = V^L^(r). It 
describes interactions between electrons with the same spin, (b) The effective Coulomb potential V^^ (r) for 
= 0, 1, 2, 3 from top to bottom. It describes interactions between electrons with the different spins. 



for M > in terms of associated Laguerre polynomials. 

A comment is in order on the form factors. It is easy to see 

GN{q)^e-"^'Bm"FN(,q). (3.31) 

Thus, Gn (q) is exponentially smaller than {q)- For instance, F|_ (g) = Gn (q) represents the 
transfer of the up-spin electron {a =t) from the K' point (r = — ) to the K point (r = +). Hence, 
such a mixing between the K and K' points is actually negligible. 

The projected Coulomb Hamiltonian is constructed by substituting the projected density (|3.27p 
into the Hamiltonian (I3.24p . 

Hn =vr j d^qV{q)pN{-Q)PN{(l) 

=''Y.Y. jd\yN-lr'xA^)D%'i-^)DiA^), (3.32) 

TT'a' XX'a' 

where ^^'^^-^-1 \x' ^^"^ effective Coulomb potential in the A^th energy level, 

yr;.,AA'(Q) = V{q)F^A-Q)Fxx'i<l) (3-33) 



with 



It is remarkable that the effective Coulomb potential depends on the spin and the valley through 
the form factors F"^,{q) characterizing Landau levels. 
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FIG. 4: (Color online) The spin dependence of the effective Coulomb potential (f), V^' {r) and V^^{r) 
in the iVth energy level for (a) iV = 0, (b) = 1, and (c) N — 2. All effective potentials agree for the 
lowest energy level {N = 0). The vertical axis is the potential in unit of e'^/AnelB- 

The typical effective Coulomb potentials in the A^th energy level at the K point are 

V^^ (q) =V{q)FN{-q)FN{q), (3.35a) 
V^^ (q) =V{q)FN-i{-q)FN-iiq), (3.35b) 
V^^ (q) =V{q)FN{-q)FN-i{q). (3.35c) 

1" I 

The potential V^^ (q) stands for the interaction between electrons with the same spin a, (q) 
for the interaction between electrons with the different spin. The spin direction is reversed at the 
K' point. We have illustrated the effective Coulomb potentials in Figsl3]and|4l 

We note that the effective potential (q) for electrons with the same spin in the A^th energy 
level has precisely the same form as in the standard QHE for electrons in the A'^th Landau level: 
Compare Figl3]Ja) with Fig. 7 in Ref.— . On the other hand the effective potential V^^ (q) for 
electrons with the different spins are entirely new. 



IV. EXCITONIC CONDENSATION 

It is necessary to pay a special attention to the zero-energy level [N = 0), since it contains both 
electrons and holes. Let us recapture the properties of the zero-energy level [Figl5]. In the absence 
of the magnetic field the band structure is given by the Dirac valleys associated with the dispersion 
relation (12. 2p . as illustrated in FigUJa). When the magnetic field is applied, the band structure is 
changed to generate Landau levels [FiglUb)]. However, the Dirac electron is subject to the intrinsic 
Zeeman effect, which induces a Landau-level mixing, and the zero-energy state emerges for up-spin 
electrons as well as down-spin holes [Figl5](c)]. There exists an attractive Coulomb force between 
an electron and a hole. Hence we expect them to make an excitonic condensation, producing a gap 
to the electron and hole states [Figl5](d)]. 

The excitonic condensation in graphene has been studied in various contexts^i^i^i^. Here, 
we present a clear-cut approach to this problem on the analogy of the BCS theory. Our physical 
picture is summarized in FiglS) The remarkable point is that the kinetic term is quenched in each 
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(a) (b) tc) m 




FIG. 5: (Color online) (a) Electron states (blue region) and hole states (red region) in the Dirac valley at the 
K point, (b) When the magnetic field is applied, the band structure is changed to generate Landau levels, 
(c) The intrinsic Zeeman effect induces a Landau-level mixing, and the zero-energy state with the up-spin 
electron and the down-spin hole emerges, (d) They make an excitonic condensation to form a BCS-type 
state, producing a gap to the electron and hole states. 



Landau level, and hence it is also absent in the zero-energy level. This simplifies considerably our 
analysis. 

"t I I 

The field operators present in the zero-energy level are 4'l_^_{x), il)l_{x) for electrons and i^l^_^_{x), 
for holes coming from the = Landau levels [FiglUa)]. It is sufficient to investigate 
electron-hole pairs at the K and K' points separately. An exciton composed of an electron and a 
hole belonging to different valleys is fragile, because their effective Coulomb potential involves the 
factor e~^^B\^\ compared with the one within the K point: See (|3.29l) and (jS.Sip . For definiteness, 
we consider the K point [Figl5]. 

The effective Coulomb potential (I3.33P is simple in the lowest Landau level, 

yeff(g) = V{q)e-^BiV\ (4.36) 
We consider the effective Hamiltonian together with this effective potential, 

H'''' = - j d^xd^V'^^ix - y)^l\x)i;l{x)^pi\y)^|;iiy), (4.37) 

where i^l{x) and V'h(^) are the up-spin electron field and the down-spin hole field, respectively. 
The Hamiltonian is rewritten as 

H''' =-^j V^^{q)^l:l\k + q)^l:\^{k')^l:l{k' - q)i;l{k). (4.38) 

Here and here after, under the symbol /, the integration over momentum variables is understood. 

We derive the gap equation, following the analysis familiar in the BCS theory. We take the 
terms satisfying q = k' — k as the dominant ones, and approximate the Hamiltonian as 

H""^ ^-Tij F^ff(fc' - fc)^Jt(fc')^it(fc')^i(/,)^T(fc). (4.39) 
We define the singlet excitonic gap function by 

A (fe) = / d^k' V'^ik' - k){i>mk')i;l\k')), (4.40) 
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which can be taken to be positive without loss of generahty. 
The mean-field Hamiltonian reads 



By setting 



ii-^ff ~ -vr / [A (k) (fc) (fc) + V'l^ (fc) < (k) A* (fc)]. 



^1 (fc) =— [Vl + 

^2 (k) =^[^r (fc)-V'i(fc)]2, 



it is easy to diagonalize (|4.41l) as 

H'^^TT j cfkA{k) [^{ (k) ^1 (fc) + (fc) (A;)]. 
Since and ^'2 are free fields, the ground state is given by solving 

^l(fc)|«'exc) = ^-2 (fc)|«5exc) =0, 



or 



\^^^^) = lll=[l + ^l^{k)^l^ (fc)]|0). 



This is a BCS-type state representing the condensation of electron-hole pairs. 

Due to the Fermi statistics the thermodynamical average (^J (k) (k)) is given by 

{^l {k)^^{k))={^{k) ^2 {k)) 



1 + e^(*=)/'=BT ' 
{k)^i{k))={^ik) M/2 (fc))=0, 
where fee is the Boltzmann factor. Combining these with (I4.42p we obtain 



i^l^ (fc)^;t(fc)) = -tanh 



A(fc) 



Substituting this into (I4.40p we derive the gap equation, 



A (fc) = i j (fk'V^^ik' - k) tanh 



A(fc') 



2kBT 

In the limit T ^ 0, the zero-momentum gap A (fc) is given by the dispersionless relation 

- Ao. 



A(k) 



T=0 



'vr e 
VTi / — - 



2 AireiB 

For finite temperature, assuming the gap is dispersionless, we obtain the relation 



At 



tanh ■ 



A7 



2A;bT 



(4.41) 

(4.42a) 
(4.42b) 

(4.43) 
(4.44) 
(4.45) 



(4.46) 



(4.47) 



(4.48) 



(4.49) 



(4.50) 



The critical temperature Tq at which At = is solved as Tq = Ao/2/cb. 

According to the diagonalized Hamiltonian (|4.43|) . the excitonic condensation provides them 
with the mass Aq, resolving the electron-hole degeneracy in the zero-energy state. Combining the 
results at the K and K' points, the 4-fold degenerate levels split into two 2-fold degenerate levels 
[FiglUJd)] with the gap energy (|4.49p . This leads to a new plateau at u = 0, as illustrated in FiglH] 
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(a)QH 

Plateau 




(b) Energy Level 



K K' 



■v=10 
■v=6 

v=2 



(c) Energy Level 



^v=0 
::iizi^v=-1 

■v=-2 

:::::: ^v= -4 

■v=-6 

v=.io^:::::::-:?. 



FIG. 6: (Color online) (a) The QH conductivity in graphene. The dotted black curve shows the sequence 
ly ~ ±2, ±6, ±10, • • • , while the soHd red curve the sequence i/ = 0, ±1, ±2, ±4, ±6, ±8, • • • . (b) The energy 
level within noninteracting theory, indicated by the number in circle. The spin is indicated by a solid red 
(open blue) arrow for electron (hole) at the K and K' points. The number attached to a solid red (open 
blue) arrow shows from which Landau level the electron (hole) comes from, (c) When Coulomb interactions 
are included, the zero-energy level splits into four nondegenerate subbands, while each nonzero-energy level 
splits into two 2-fold degenerate subbands with the exact U(l) symmetry. 



V. VALLEY POLARIZATION 

It is now possible to treat electrons and holes separately since the gap has opened between the 
electron and hole bands. Furthermore, it is enough to study only electrons due to the electron-hole 
symmetry. We show that the Coulomb effect modifies the energy spectrum of the noninteracting 
theory so that plateaux emerges at u = ±1, ±2n with n = 1, 2, 3, • • • [FiglB]. 

We express the Coulomb Hamiltonian (I3.32P as 

Hn = 5^y;;;:^;^:'"'<t(m)<,(n)cf (,)c^;(j), (5.1) 
where the summations over repeated indices, m, n, z, j, o", a' , r, r'. A, A' are understood, and 
y^j"^ ~ I d'qV{q)F^^,i-q)F^,y{q) 

x(m|e-*[-''+-^-^'^l^|n)(z|e-['^+"^-^'^]^|j). (5.2) 

We introduce 

T/rr'AA' _ t .TT;Tr'AA' -r.rr'AA' _ t .TT;Tr'AA' n\ 
VN;B - ^N;iijj ' ^N;X - ^N;ijji ■ \^-^} 

ij ij 

They represent the direct and exchange Coulomb energies, respectively. It is easy to evaluate them 
numerically: See FiglT^a) for V^^^ . We analyze the lowest energy level and higher energy levels 
separately. 
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FIG. 7: (a) The exchange energy of the TVth energy level as a function of N . The vertical axis is 

the energy in unit of e^/47re£_B. Note that V^.x~" = '^N-vy^- (b) The function tan^aAr given by (|5.12p . 
It approaches 1 assymptotically, liniAr^oo tan^ ofAr = 1. 

A. Lowest Energy Level 

We first study the lowest energy level, which contain two degenerate states [FiglUb)]. The 
degeneracy is resolved obviously by a strong extrinsic Zeeman effect if it exists. We consider the 
case where the extrinsic Zeeman energy is absent at temperature T = 0, or is present but quite 
small compared with the thermal energy at T / 0. We ask whether the degeneracy is resolved even 
in such cases. 

As a generic trial function we take 

1%) = n (^c^l" (n) + vc^l (n)) |0) , (5.4) 

n 

with |np + = 1. The state covers the entire SU(2) space of the lowest energy level for electrons, 
when the two parameters u and v are varied. 

We calculate the Coulomb energy (i^)o = (^o|-f^|^o) with (|5.ip . Most terms are SU(2) invariant, 
but there exists a noninvariant term, 

(^fnonmv>0 = 27r\uv\\++++ . (5.5) 

Hence the SU(2) symmetry is broken explicitly into the U(l) symmetry by Coulomb interactions. 
The energy is minimized either u = oi v = due to the term (|5.5p . corresponding to the state 

n n 

The ground state is either |$q) or though there exists still the Z2 symmetry. The energy 

barrier is of the order of the Coulomb energy, which is much larger than the thermal energy. 

We conclude as follows: The two levels split explicitly by an extrinsic Zeeman effect if exits. 
Then the first energy level is up-spin polarized, and the second energy level is down-spin polarized. 
Even without such an extrinsic Zeeman effect, driven by the Coulomb exchange interaction, the 
spontaneous breakdown of the Z2 symmetry turns the system into a QH ferromagne t^^'^^ . It is 
reasonable to call it the Ising QH ferromagnet due to the Z2 symmetry. In any case, a plateau 
emerges at u = 1, where the activation energy is of the order of the typical Coulomb energy as in 
the conventional QHE2^>2Z. 
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B. Nth Energy Level 



We next study the A^th energy level with > 1. It contains four degenerate states with the 
SU(4) symmetry in noninteracting theory. However, the projected density is invariant only under 
U(l)xU(l)xZ2. Hence we take a set of trial functions by requiring this symmetry, 



= n (sinoTve^^ cj] (n) + cosQ^e"*^ (n)j |0) , (5.7a) 

n 

l^iv) = n (cosQ^e~*^-^ci! (n) + sina^e^^^ci^ (n)) |0) , (5.7b) 

n 

where the phase factors e*^^ and e*^^ assure the U(l)xU(l) symmetry. These two states are degen- 
erate, ($]y|i?iv|$jY) = (^>jy|iJAr|$jy), due to the Z2 symmetry. 

It is easy to determine the angle ajy by minimizing the Coulomb energy (<I>^|i/iv|$^). After 
some calculations we find 



{^^jylHNl^'^jy) = A\ sinQTvl^ + B\ cosaAr|^ + C| sinQArcosaAr|^, (5.^ 



where 



A =^^;i+^ - ^A^;x^^' B = y^--- - y^---, (5.9a) 

C =V^:^- + V^-^-^+ - - (5.9b) 

with ([531) • Here, Vj^^'^^ = V^.^ = V^.+ = ^^7.^++ since F^^,{0) = 1. Furthermore, because 
of (I3.3ip . and V^-x^ are exponentially smaller in {a/H.B)'^ than V^^^ or V^.x , and 

can be neglected. Hence 

2 B-CI2 Vn^^' 

Now it follows from that V^-x"" = V^^i^x • We can also see [FigllKa)] 

implying that the Coulomb energy of an electron in higher Landau level is lower. It follows that 

Y++++ 

tan^aAf ^ ++++ > 1- (5.12) 

We have depicted tan^ as a function of N in FiglUb). 

We conclude that, if we take the state l^jy), more electrons are present in the Dirac valley at the 
K point than at the K' point since | sinaAr| > | cosaAr|. Namely, the valley polarization has occurred 
both in and \<^\^). This can be understood physically as follows [FigHlJa)]: Up-spin electrons 
in the K point belong to the (A^-l)th Landau level but those in the K' point belong to the A^th 
Landau level. It is easier to fill Landau sites at the K point because the Coulomb energy is lower 
in higher Landau levels. The valley polarization disappears as A^ — > 00, since V^^{7x = ^^v x 
the limit [FigllKb)]. 

An extrinsic Zeeman effect open a gap between these two spin polarized states. Even without 
such an effect, driven by the Coulomb exchange interaction, the spontaneous breakdown of the Z2 
symmetry turns the system into a QH ferromagnet. Note that it has still the 2-fold degeneracy. In 
any case, a plateau emerges at u = AN [FigEfc)], where the activation energy is of the order of the 
typical Coulomb energy as in the conventional QHE2^>2I. 
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VI. MULTILAYER GRAPHENE SYSTEMS 



A. Bilayer Graphene (Bernal Stacking) 

We proceed to generalize the above analysis to a bilayer graphene, which is a system made of 
two coupled hexagonal lattices according to the Bernal stacking. In the absence of magnetic field, 
the low-energy spectrum of the bilayer graphene is known^i^>^i^ to be parabolic, 

£{k)(x\k\'^, (6.1) 

near the K and K' points. In the presence of the magnetic field, we can reformulate the model 
Hamiltoniani^>^>^ as the generalized Dirac Hamiltonian defined by 

= diag. - ^/QlQ±j , (6.2) 

together with 





Here, A = hwcO?, with a given by (|2.14l) . We also consider the generalized Pauli Hamiltonian 

H+ =Q+Q+ = {huJcf I I , (6.4a) 

H- =Q-Q- = {huj.f I I . (6.4b) 



We may switch off the magnetic field in this formula, and reproduce the energy spectrum (|6.1jl . 
The eigenvalue of the Hamiltonian is derived as 

H^\N) = [£^\£p,-£^\-£p) \N) (6.5) 

with the eigenstate |A^) = {N\)-'^/^{a^)^\0), where 

£p = £^^ = nwcVN{N-l), (6.6a) 

£+^ = £-^ = nwc^/{N + 2) {N + l) (6.6b) 

for = 0, 1, 2, 3, as illustrated in FiglEl It is interesting that two Landau levels mix to create one 
nonzero-energy level, but that four Landau levels mix to create the zero-energy level. Thus there 
exists the 4-fold degeneracy in the nonzero-energy level but the 8-fold degeneracy in the zero-energy 
state, as results in the bold-face series (11.21) . This agrees with the previous result^i^i^. 

We include Coulomb interactions. The A^th energy level {N / 0) has the same structure as 
in the monolayer graphene system. Coulomb interactions make each energy level split into two 
subbands. 
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FIG. 8: (Color online) (a) The QH conductivity in bilayer graphene. The dotted black curve shows the 
sequence ly = ±4, ±8, ±12, • • • , while the soHd red curve the sequence ly = 0, ±1, ±2, ±3, ±4, ±6, ±8, • • • . 
(b) The energy level within noninteracting theory. The spin is indicated by a solid red (open blue) arrow 
for electron (hole) at the K and K' points. The number attached to a solid red (open blue) arrow shows 
from which Landau level the electron (hole) comes from, (c) When Coulomb interactions are included, the 
zero-energy level splits into eight nondegenerate subbands, while each nonzero-energy level splits into two 
2-fold degenerate subbands with the exact U(l) symmetry. 



We discuss the zero-energy state in some details. There are electron-hole pairs coming from the 
= Landau level and the = 1 Landau level for electrons and holes [Figl8](b)]. At the K point 
there are four gap functions Aqo (fc), An (fc), Aqi (k) and Aio (fc), 



(fc) = / d'k'V^^ik' - fc)(V'fTt(fc')^A^'it(fc')), 



(6.7) 



where •0^^^ and ■0^ ' ' creation operators of electrons in the A^th Landau level and holes in the 
A'^'th Landau level, respectively. Repeating similar analysis we have made in Section ITVl we solve 
the gap equations at T = as 



Aoo (fc) =Ao, 

An (fc) =^Ao, 

Aoi (fc) =Aio (fe) = ^Ao, 



(6.8a) 
(6.8b) 
(6.8c) 



where Aq is given by (14491) : Aq = iry^TrJlie^ /AneeB). Thus, the 8-fold degenerate zero-energy 
level splits into four 2-fold degenerated subbands, producing plateaux at = ±2, ±4. Furthermore 
we have Ising QH ferromagnets at v = ±l,ib3. Consequently, we predict the full series (11.21) . 



B. Trilayer Graphene (Rhombohedral Stacking) 

The above analysis is applicable also to a trilayer graphene with the ABC stacking (rhombohedral 
stacking). In the absence of the magnetic field, the low-energy spectrum of trilayer graphene has 
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FIG. 9: (Color online) (a) The QH conductivity in trilayer graphene. The dotted black curve shows the 
sequence v = ±6, ±10, ±14, • • • , while the solid red curve the sequence z/ = 0, ±1, ±2, ±3, ±4, ±5, ±6, 
±8, ±10, • • • . (b) The energy level within noninteracting theory. The spin is indicated by a solid red (open 
blue) arrow for electron (hole) at the K and K' points. The number attached to a solid red (open blue) 
arrow shows from which Landau level the electron (hole) comes from, (c) When Coulomb interactions are 
included, the zero-energy level splits into twelve nondegenerate subbands, while each nonzero-energy level 
splits into two 2-fold degenerate subbands with the exact U(l) symmetry. 



been argued^^ to be cubic, 

£{k)oc\k\^, (6.9) 

near the K and K' points. In the presence of the magnetic field, we can reformulate the model 
Hamiltonian^^ as the generalized Dirac Hamiltonian as 

= diag. (VOlgl, -x/Qigl) , (6.10) 

together with 



Q+=\ , Q- = I . (6.11) 





Here, A = hwcO? , where a given by (I2.14p . 

The eigenvalue of the Hamiltonian is derived as in (16. 5p with 

= ^N^ = nucVNiN -l)iN -2), (6.12a) 

^N^ = = f^cViN + ^) iN + 2) (iV + 1) (6.12b) 

for = 0, 1,2,3, as illustrated in FiglS) There exists the 4-fold degeneracy in the nonzero-energy 
level but the 12-fold degeneracy in the zero-energy state, as results in the bold- face series (|1.3p . 
This agrees with the previous resulti^. 

We briefly argue how the energy spectrum is modified by Coulomb interactions. The effective 
Coulomb potential depends on the spin and the valley degree of freedom precisely by the same 
mechanism. The A^th energy level (A^ / 0) has the same structure as in the monolayer graphene 
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system. Hence, each energy level splits into two subbands. The ground state is describe by a 
formula similar to (15. 7|) . where the valley polarization is occurred. There are electron-hole pairs 
in the zero-energy state, which make excitonic condensation by the same mechanism as in the 
monolayer and bilayer cases. Thus, the 12-fold degenerate zero-energy level splits into six 2-fold 
degenerated subbands, producing plateaux at u = ±2, ±4, ±6. Furthermore we have Ising QH 
ferromagnets at = ±l,ib3, ±5. Consequently, we predict the full series (I1.3I1 . 

VII. DISCUSSIONS 

The most intriguing property of the graphene system is that the intrinsic Zeeman energy is 
precisely one half of the cyclotron energy for electrons and holes. It leads to the symmetry group 
SU(4) in the noninteracting theory, where the Hall plateau emerges at = ±2, ±6, ±10, • • • . This 
series is the first experimental resulii^ of the QHE in graphene. When Coulomb interactions are 
included, the symmetry SU(4) is broken so that the Hall plateau emerges at u = 0, ±1, ±4, ±8, • • • . 
This series has been found experimentally^ when larger magnetic field is applied. We have shown 
that the Ising QH ferromagnets appear at v = ±1 due to a BCS-type condensation of electron-hole 
pairs. 

We have emphasized that one energy level contains up-spin and down-spin electrons coming from 
two neighboring Landau levels. Since Coulomb interactions are different for electrons in different 
Landau levels, we have derived a remarkable consequence that the effective Coulomb potential 
depends on the spin and the valley degree of freedom [FigslS] and |4]. As a result, the valley 
polarization is occurred on the ground state. We wish to explore new physics associated with this 
peculiar Coulomb interaction in forthcoming papers. 
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APPENDIX A: DIRAC ELECTRONS IN GRAPHENE 

In this appendix, we derive the second-quantized Dirac Hamiltonian (12. 3p together with (|2.4p 
for electrons in graphene. We start with a study on the quantum-mechanical states of electrons 
based on the k • p approximation^i^^. 

First of all, the wave functions fg{x) are given in terms of envelope functions Fg{x) as^ 



where S is the site index {S =A,B) and r is the valley index (r =K,K'). The envelope functions are 
normalized as \F^{x)\'^ -\- \F^'{x)\'^ = 1/2 for S =A,B. They satisfy the Schrodinger equation^! 



ff{x) = e'^'-^Ffix) 




(1.13) 



HoF{x) = £F{x) 



(1.14) 
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where S is the eigenvalue, and the Hamiltonian Hq is the 4x4 matrix operator, 



Px - iPy 
Px + iPy 



v 


















Px + iPy 



with pk 



Px - iPy 

-ihdk- Since this is block diagonal, it is convenient to set 



(1.15) 



= vpia- -p) = vf 




(1.16a) 



(1.16b) 



with p' = {px, —Py)- Note that = axH^ax, where ax is the generator of the mirror symmetry. 
The eigenfunctions of (r =K,K') are given by^i 



(1.17) 



with the eigenvalue £{k') = ahvp\k'\, where a stands for the helicity, a = ±. Note that hk' is the 
momentum of electrons, \k'\ <C \K\ measured from the K or K' point. 

Let us explain the notations we have used in (|1.16p and p.l7|l . We recall that the cr • p is 



the helicity operator except for the positive normalization factor. Thus, F^'^{x) has the positive 
helicity and a positive energy, while F^'^{x) has the negative helicity and a negative energy. Since 
the energy and the helicity have the same sign for the envelope function -F^'^(a;), it describes the 
right-handed Weyl fermion by definition. Hence we have put the index "R". Similarly we have 
assigned the right-handed Weyl fermion at the K' point. 

There are only four independent two-component envelope functions given by (I1.17P with r =K,K' 
and a = ±. However, there are more quantum-mechanical states for electrons. They are the 
chiral symmetric copies; the left-handed Weyl fermions F^''^{x) at the K point and F^'^{x) at 
the K' point. The corresponding Hamiltonians and envelope functions are constructed by the 
chiral transformation generated by the Pauli matrix o"^ as = azH^az = —VF{a -p), H^' = 
^'^ — —VF{a • p'), and 



X 



:it:R, 



X) 




(1.18) 



Note the the energy and the helicity have the opposite sign for the left-handed Weyl fermion. 

The energy spectrum is symmetric between the positive and negative energy states. There exists 
one electron per one carbon and the band-filling factor is 1/2 in graphene. Namely, all negative- 
energy states are filled up, as is a reminiscence of the Dirac sea. Hence, we have electrons and holes 
as physical excitations. 

In this way there are eight types of quantum-mechanical states for electrons in graphene, cor- 
responding to the spin degree of freedom (cj = it), the electron-hole degree of freedom and the 
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valley degree of freedom (r =K,K'). To carry out the second quantization it is necessary to use 
all these quantum-mechanical states together with the four Hamiltonians. We arrange them as 
= vpa • P75, H^' = —vpa • p'75 where we have introduced the Dirac 75 matrix in the Weyl 
representation, 




(1.19) 



They are summarized into the Dirac Hamiltonian (|2.4p in text, or 

i^D = Vp {raxPx + (^yPy) 75> (1-20) 

where r = + (— ) for the K (K') point. 

According to the standard prescription the second-quantized Hamiltonian is constructed as in 
(|2.3p . 01 H = X]r=± / '^t{x)H^"^T-{x) , where is the quantum mechanical Hamiltonian. The 
field operator is expanded as "^rix) = ^'er(a^) + ^hrl^j), with 

*h.(^) = E / ^dp{k')vl{k')e-^''''^. (1.21b) 

We have introduced the annihilation operator c"{k') of an electron with the eigenfunction u'^{k'), 
and the creation operator dr^k') of a hole with the eigenfunction v'^{k'), where we have set u'!^{k') = 
F^_^^{k'), vL{k') = FL'^(fc'), vl{k') = F^'^(-fc') and wL(fc') = FL'^(-fc') in accord with the 
standard notation in the Dirac theory. In passing, the field operator V'r(^) fo^' electrons in graphene 
is to be constructed with the use of the wavefuncion (|1.13p . and hence is given by (|2.7p in text. 



II. LANDAU-LEVEL PROJECTION 



In this appendix we derive the formula (I3.27P for the projected density operator Pn{q)- We 
first review how the Landau-level projection is made in the conventional QHE, where the density 
operator is given by p{x) = ip^ {x)%l^{x). We consider electrons confined to the A^th Landau level, 
where Fock states are given by (I2.22p . The field operator is expanded as 

tP^x) = '^{x\N,n)c{n), (2.22) 

n 

with c(n) the annihilation operator of electrons acting on the state \N,n), {c{n),c^{m)} = 6nm- 
The projected density operator is 

Pj<^{x) =Tp\x)i;{x) = ^{N,m\x){x\N,n)c\m)c{n). (2.23) 

mn 

Its Fourier transformation is 

p^{q) = —y" [(fxe-''i''{N,m\x){x\N,n)c^m)c{n). (2.24) 

mn 
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Here, we decompose the coordinate into the guiding center and the relative coordinate, x = X + R, 
where X and R act on the Fock states \n) and |A^), respectively. Since they commute each other, 
we obtain 

p^{q) = -1 y{N\e-'''^\N){m\e-'''^\n)c^m)c{n), (2.25) 

mn 

or Pi^iq) = FN{q)p{q), with 

FN{q) = (iV|e-*'?^|iV), (2.26a) 

p{q) = ^ V(m|e"^''^|n)ct(m)c(n). (2.26b) 
27r ^ — ^ 

mn 

We call Fp^{q) the Landau-level form factor, and p{q) the bare density operator. 

In the graphene QHE one energy level contain four different types of electrons described by 
(|2.23ll . For instance, the projected density operator for up-spin electrons is given by 

p\{x) ={N, m\x) {x\N, n)c^_^ {m)c^_^{n) + {N — l,m\x) {x\N — 1, n)cl^(m)cl(n) 

+ e-'^'^''{N,m\x){x\N - l,n)c^^{m)cl{n) + e^'^''{N - l,m\x){x\N,n)cl! {m)cl{n), 

where the summation over the indices n and m is understood. By repeating the above process, this 
reads 

pUq) =Fn% (9) i^lV (q) + Fn-i,N-i (9) Dll (q) 

+ Fn^n-i (9) (g) + Fn-i,n (g) ^-H- (q) > (2-27) 

where 

F]7'm (q) = (iV|e-*['?+-^--'^]«|M), (2.28) 

and 

D^f,{q) = ^ y (m|e-^[''+"-^-"'^^|n)<t(^)c-;(,„). (2.29) 

mn 

A similarly formula is derived also for p\^{q). Adding p\{q) and p\q{q) we obtain (I3.27P in text. 
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